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Abstract 

In this article we consider a simple model for a c-axis Pb/YBa2Cu307_(5 
Josephson junction. The observation of a nonzero current in such a junc- 
tion by Sun et al. [A. G. Sun, D. A. Gajewski, M. B. Maple, R. C. Dynes, 
Phys. Rev. Lett. 72, 2267 (1994)] has been taken as evidence against d-wave 
superconductivity in YBa2Cu307_5. We suggest, however, that the pairing 
interaction in the Cu02 planes may well be d~wave but that the CuO chains 
destroy the tetragonal symmetry of the system. We examine two ways in 
which this happens. In a simple model of an incoherent junction, the chains 
distort the superconducting condensate away from d^2_y2 symmetry. In a 
specular junction the chains destroy the tetragonal symmetry of the tunnel- 
ing matrix element. In either case, the loss of tetragonal symmetry results 
in a finite Josephson current. Our calculated values of the critical current 
for specular junctions are in good agreement with the results of Sun and 

co-workers. 
74.50.-hr, 74.20.Fg 
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I. INTRODUCTION 



The debate over the symmetry of the order parameter in the high copper-oxide su- 
perconductors has intensified over the last few years because of a number of suggestive ex- 
perimental findings-lll The discovery of linear low temperature behaviour in the penetration 
depths of single crystals of YBa2Cu307_5 (YBCO)i'i and Tl2CaBa2Cu208„5@ has been taken 
as support for rf-wave super conduct ivityil, although such experiments have been unable to 
exclude anisotropic s-wave models.i^ NMR relaxation ratesi0 have been interpreted in 
terms of a (i-wave order parameter.0lli More recently, angle resolved photoemission ex- 
periments (ARPES)EHi have been able to map out the Fermi surface in the normal and 
superconducting states for Bi2Sr2CaCu208+5 and have found an anisotropic gap with nodes 
(or at least minima smaller than the resolution of the experiment) located approximately 
along the diagonals of the Brillouin zone. This is strongly suggestive of a rf-wave gap al- 
though the experiments have been criticized because they sample only the surface states of 
the crystals. 

A recent generation of experiments has attempted to resolve the issue of order param- 
eter symmetry by measuring its relative phase between different regions of the Brillouin 
zone. These experiments are all based on the fact that the current through a Josephson 
junction depends on the phase difference between the condensates on either side. Some 
of the experiments which have been performedEUll have attempted to measure the phase 
difference between electrons tunneling through different faces of a single crystal of YBCO. 
In a (i-wave material, these electrons will have a phase of vr relative to each other. Two 
of the experiment suggest that YBCO has a rf-wave order parameter while the thirdii 
suggests an s-wave order parameter. More recently, there have been experiments which 
attempt to detect phase shifts of vr across YBCO/YBCO junction^li3 in which the crystals 
are misaligned. These phase shifts, which are a signature of rf-wave symmetry, have been 
found in both of the cited experiments. 

The experiment in which we are interested in this article is that of Sun et a/.0. It is 
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slightly different than the others: it relies on the fact that in a c-axis Pb/YBCO junction 
(in which the junction face is perpendicular to the YBCO c-axis) the Josephson current will 
vanish if YBCO has a d-wave gap. Since the experiment finds a small, but finite Josephson 
current {JcRn = 0.3 — 0.9 mV, where -R„ is the normal state resistance of the junction and 
Jc is the critical Josephson current) the authors conclude that the order parameter cannot 
be purely c/-wave. On the other hand, their results are also inconsistent with simple s- 
wave theory since the measured critical voltages JcRn are an order of magnitude lower than 
expected. Following Ambegoakar and Baratoffi, and assuming that the gaps in the YBCO 
and Pb are Ay ~ 14meV and Ap ~ 1.4meV respectively, they find 



" " e Ay + Ap" \^ Ay + Ap / 
~ 8 mV 

where K is the complete elliptic integral of the first kind. 

In fact, the experiment of Sun et a/.0 does not immediately rule out YBCO having a 
dx2-y2 order parameter. Tanakali has shown, that while the usual treatment of the barrier 
as a second order perturbation does lead to a vanishing Josephson current, higher order 
terms in the perturbation series will not vanish. Unfortunately, attempts by one of the 
current authorJ^^ to fit the experimental results with a fourth order calculation have been 
unsuccessful. On the other hand, it may be reasonable to expect that YBCO does not have a 
gap structure with a full d^2_y2 symmetry since (optimally doped) YBCO is not tetragonal. 
Band structure calculationsHS show that the CuO chains both contribute a piece of Fermi 
surface with orthorhombic symmetry and distort the band structure in the Cu02 planes. 
O 'Donovan et a/.0 have a simple model of this: they consider a single Cu02 plane in which 
the Fermi surface is distorted slightly away from tetragonal symmetry. They find that the 
reduced symmetry is strongly reflected in the gap, which picks up a small extended s-wave 
component while retaining its nodal structure. Another point of view is that of Xu et a/.S, 
who treat YBCO as having tetragonal symmetry but take the gap to be s + id. Despite the 
very different starting points, both of these articles reach the conclusion that a gap which 



is roughly ten percent s-wave will explain the results of Sun et al. 

In this article, we suggest that the nonzero critical current seen in Pb/YBCO c-axis 
Josephson junctions may be due to the CuO chains. We consider a phenomenological model 
for YBCO in which the unit cell contains a Cu02 plane and a CuO chain and calculate 
the tunneling current for a Pb/YBCO junction as a function of the chain-plane coupling 
strength. In Sec. |l|, we introduce our model for YBCO. In section |T| we derive an expression 



for the critical Josephson current. In Sec. |I^ we discuss the results of numerical calculations. 
We finish with a brief conclusion in Sec. 0. 

II. MODEL 

We wish to consider a simple model for a Pb/YBCO Josephson junction. For the sake of 
clarity, we will take the Pb (YBCO) to be on the left (right) side of the junction. We treat 
the Pb as an ordinary free electron metal with an isotropic BCS gap. The YBCO is treated 
with a simplified model in which there are alternating layers of chains and planes. The 
planes and chains are weakly coupled by coherent electron hopping and the planes contain a 
BCS-like pairing interaction with a dx2_y2 symmetry. The chains are driven superconducting 
by a proximity effect. This model is related to models studied elsewhere&B in which both 



layers are treated as planes and (with the exception of Ref. [ |3^ where it is d-wave ) the 
order parameter is s-wave. It should be emphasized that this model is only suitable for 
weak chain-plane coupling. As the coupling is increased, the pairing potential begins to 
affect electrons in the chains directly. The problem becomes more complicated in this case 
and has only been examined in various special limits. 

The experiment we wish to describe is one in which the c-axis of the YBCO is normal to 
the junction.0 Furthermore, we assume that the tunneling junction is adjacent to a Cu02 
planeii. Then our Hamiltonian can be written: 

H = Ho' + Ho" + T (1) 

where Ho' and Ho*^ describe the uncoupled Pb and YBCO subsystems and T describes the 



coupling through the junction. We define Cko-, aiko- and a2k<T, to be the electron annihila- 
tion operators with wavevector k and spin a, in the Pb, YBCO planes and YBCO chains 
respectively. We can write 

W-NV = E<k)cLck. 

- Y: [A'4^cLki + A'*c_kiCkT] (2a) 

k 

Ho'' - N> = J2 [alkaaika6(k) + a^kaa2ka6(k) 

k,(T 

+ E [^(k)aL.a2k. + t*(k)a5k<,aik<.] (2b) 

k,(T 

- E [^kalkT^Lki + A^ai-kiaikT 

k 

T = E ^kqCko-aiqcr + Tj^qajq^Ckcr (2c) 
k,q 

The normal state dispersion in the Pb is e = h^k"^ /2m* — /i, and the mean field order 
parameter, A', is isotropic. The dispersions and ^2 are the normal state dispersions of 
the YBCO planes and chains in the limit of no chain-plane coupling: 

(^1 = — 2ai[cos(A;^) + cos(A;j^)] — /ii (3a) 

^2 = -2(72 cos{ky) - /i2 (3b) 

with — TT < kx,ky < TT and where fii and fi2 include information about both the chemical 
potential and the offset of the bands from one another. The strength of the plane-chain 
coupling is given by t(k) which, in the tight binding limit, dependsil only on kz- To simplify 
matters further, we take the chain-plane distances to be the same on either side of a chain, 
so that 

t{kz) =tocos{kj2), (4) 

with — TT < kz < IT. The mean field order parameter in the YBCO, A^, is a thermal average 
of electron pairs in the CUO2 planes only. This is because we make the ansatz that the 
pairing is localised to the planes: 



(5) 



where Vkk' is the BCS-hke pairing interaction in the planes. If, for simphcity, we assume 
that the pairing interaction is separable, so that Vkk' — VVkVii'^ with r/k = cos(A;3;) — cos(A;y) 
for a d-wave interaction, then we may write AJ^ = AgT^k- 

We will find it convenient to work within the Nambu formalism, in which we define 



C(k) 



and 



A(k) 



CkT 



aikT 



^l-ki 



a-2kT 



*2-ki 



so that the uncoupled Hamiltonians may be written. 



H« = ^Ct(k)7i|,(k)C(k) + i:At(k)7iS(k)A(k). 



(6a) 



(6b) 



(7a) 



with 



and 



^o(k) 



e(k) -A' 
-A'* e(-k) 



njr 
/to 



6(k) -Ai t{k) 

-A[* -6(-k) -r(-k) 
r(k) e2(k) 
_i(_k) -6(-k) 

The eigenvalues of the Hamiltonian matrices are El^^{k) with El — —E'-), 



E^ = V^2^^' 
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(7b) 



(7c) 



(8a) 



and El = (£;;, EL, -EL, -E^), 



E'^ 



\ 



+ ^^[(ei + e2)^ + Ai;^]. 



and the unitary transformations which diagonahse the Hamiltonian are 



^ndUl^=Ul{E^), 



W'(k) 



tA 
tB 



B^e-{A^C + E^^-i,){E^^-i^ 



The Hamiltonian may, therefore, be written 

Ho-N/. = ^^C|(k)Q(k)£;^(k) 

k i=l 

+ EEAl(k)A,(k)£;[(k) 



k i=l 



where, for example, Ci(k) = X^j Z4j(k)Cj(k). 

In the normal state, the YBCO band energies are 



. iLlia ± /(O^. (u) 

and their Fermi surfaces are given by = .^i^2- In Fig- we plot the Fermi surfaces in the 
kxky plane for a range of t. When t = 0, the two pieces of Fermi surface are just the Fermi 
surfaces of the isolated planes and chains, given by = and ^2 = respectively. There 
is a Fermi surface crossing at = ^2 = 0. As t is increased, the Fermi surfaces are pushed 
apart, and the crossing becomes an avoided crossing. Far away from the avoided crossing, 
the each piece of Fermi surface is predominantly chain or plane in character. However, near 
the avoided crossing, the two bands are hybridizations of the chain and plane states. In 
the superconducting state this has the important effect of distorting the gap away from the 
dx2_y2 symmetry of the pairing interaction. In Figs, ^(b) and (c), we show the quasiparticle 
band structure [given by Eq. in the superconducting state along the two paths ky = 
and ky = Skr^/S (with kz = in both cases). For comparison purposes, we show the same 
spectra in the t = limit. Along ky = 0, the two pieces of Fermi surface are far enough apart 
that one piece is predominantly chain-like while the second is predominantly plane-like. If 
we take the term "gap" to mean a local minimum in along paths of the type ky = ak^ 
then we can see that there is a double gap structure along the ky = direction. The larger 
of the two gaps can be identified with the CUO2 plane and is perturbed from |Ak| by a term 
of order | Ak|t^/ {C,i — ^2)- The second gap is the induced gap in the chains, and it is of order 
|Ak|tV(Ci -"^2)- The second path, ky = 8k^/3, passes through the avoided crossing. In Fig. 
^( c ), the band structure near the Fermi surface has little in common with the band structure 
in the t = limit, and the two gaps are nearly equal to each other, but very different from 
Ak. It is the effect of the chain-plane coupling in this region of the Brillouin zone which 
produces the finite c-axis Josephson current. 

The shape of the Fermi surfaces in Fig. ^(a) was chosen to qualitatively resemble the 
results of first principles band structure calculations. 00 Such calculations find a Fermi 
surface that has four pieces, two of which are similar to the ones shown in Fig. ^(a). The 
remaining two pieces of Fermi surface have a (nearly) tetragonal symmetry. These have not 
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been accounted for here since the goal is to describe the effects of orthorhombic distortion 
with a simple model. 



III. JOSEPHSON CURRENT 



The current generated by T is eN', where N' = I^k.o- cj^^Cko- is the number of electrons 
in the Pb and N' = -i/h[N'-, T] so that 



(N') = 2Im Y: T^Ma^i^^)- 



(12) 



k,q,(T 



Taking T as a perturbation we find that, to lowest order, the Josephson current i 



-2e 



i_ _ J—oo 



(13) 



k,q,(T 



where rj is positive and vanishingly small, and the expectation value is now taken with 
respect to the uncoupled system. Equation (0) has both a supercurrent contribution, which 
depends on the expectation values {cl^■^(t)c^_^^^(t')) and (ai_ki(^')^ikT('^))5 and a quasiparticle 
contribution, which depends on (cj^^(t)cko-(t')) and (a|j^^(t)aiko-(t')). In our case, the voltage 
across the junction is zero and the quasiparticle part vanishes. The supercurrent can be 
evaluated by rewriting the electron creation and annihilation operators in terms of the 
superconducting quasiparticle operators and noting that, for example. 



C^(k,t) = e-^™«*/^C„,(k,0). 



(14) 



It follows directly that 

J = |lm E T,,r,,_,^|:(k)^^.(k)^{,(q)ZY-(q) ^^^J^j|j (15) 

This is the basic equation for the Josephson current in the absence of an external voltage. 
For our particular model, this expression becomes 



4e, 



im[A'* Ay r - ^. 



El 



1 



Elf{E') - E^fiE": 



E^ + E\ 



+ 2 



- E[ 



r2 



EL 



1 

+ E': 



EU{E')-E'f{EL) 

02 _ 



(16) 



In this equation k is the variable of integration for all terms associated with the Pb {ie. all 
variables with superscript I) and q is the variable of integration for all terms associated with 
the YBCO. We have also used Tkq = T*]^_^, which follows from time reversal symmetry. 
The phase is the complex phase of A'*Aq. Equation (|1^) can be written J{<f)) = Jc sin(0), 
which defines the critical current Jc- 

Although Eq. ( |T6|) is complicated in appearance, its behaviour can actually be understood 
fairly easily. First of all, since > \C,2 \ everywhere (see, eg., Fig. |l|) the sign of the coefficient 
of the first term is the same as the sign of Aq. Furthermore, the term inside the square 
brackets can easily be seen to be positive and decreasing with increasing T. A similar 
argument holds for the second term in Eq. (|16]), except that E"^ — ^| changes sign between 
different regions of the Brillouin zone. In Fig. |l|, however, we can see that E'^ — ^| < near 
the Fermi surface, so that at low enough temperatures (recall that T is less than one tenth 
of the Tc of YBCO here) the sign of the temperature dependent part of the second term 
is also determined by Aq. Whether is an increasing or decreasing function of T, then, 
depends on the strength with which the integrand contributes to the integral in different 
regions of the Brillouin zone. 

It is simple to show that, in the limit ^ 0, Eq. ([TBp becomes the well known equation 
of Ambegoakar and Baratofiii 

J-J-1.\Tk^\ Ei^r (17) 

k,q 

^ \ f{E^)-f{E^) f{E^)-f{-En \ 
\ E^ -E"- + E"- J ' 

with E^ = [^l + Aq^]-"^/^. Provided that |Tkq| is invariant under n/2 rotations in the QxQy 

plane, the Josephson current will vanish for a (i-wave order parameter. 

In order to proceed with Eq. (0) it is necessary to choose a form for the tunneling matrix 

element. The two common choices are 

l^k.qP = (18a) 
(which describes an incoherent tunneling process) and 
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|Tk,qP = -^kUkK(q)|4,„,|, (18b) 

(which describes a specular tunneUng proces^il). In Eq. ( |18b| ), P is the probabihty of 
transmission through the barrier for a single electron, and U are the thicknesses of the 
Pb and YBCO perpendicular to the junction and and are the semiclassical electron 
velocities in the 2;-direction: = de/dkz and 

de± 



I T I 



dqz 



\to sin(g,/2) | , (19) 
^(^1 - 6)^ + 4t(g.)2 



The (5-function in Eq. (|18b| ) conserves the momentum parallel to the junction face. 

In the case of specular tunneling, the choice of vl plays an important role in determining 
the magnitude of the Josephson current. For a single band material, in which there are only 
Cu02 planes, = tosm{q;^) so that the entire Fermi surface contributes to the tunneling 
process, and the Josephson current vanishes because of the antisymmetry of the d-wave 
order parameter. In Eq. (0), however, there is a weighting factor which is only appreciable 
in regions where |^i — ^2! < 2|t|. Physically, this means that currents can only flow along 
the 2;-axis in regions of the Brillouin zone near to where the Fermi surfaces cross: electrons 
travelling in the z-direction must hop between the chains and planes and, since the chain- 
plane coupling is coherent (conserves q), hopping can only take place in regions where the 
chain and plane Fermi surfaces are close together. In the q-space integral in Eq. ([T6|), 
has the effect of restricting the integral to one small region of the Brillouin zone over which 
the order parameter is roughly constant. Because of this Eq. ([T6|) is not able to distinguish 
whether YBCO is s-wave or d-wave for small to- 

Another useful way of looking at vl is that it destroys the antisymmetry of the integrand 
under rotations of tt/2. For small to, the Josephson current is approximately given by Eq. 
(p!7|). The integral is non-vanishing, however, because of |Tkq|'s lack of symmetry. 

We will finish this section with a brief derivation of the normal state junction reistance 
Rn, which is necessary to determine the critical voltage, JcRn- The calculation is similar to 
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the one performed above for the supercurrent and, as before, we begin with Eq. (|T3D. In this 
case, however, we are finding the quasiparticle current, J„, driven through the junction in 
the normal state by a vohage V, and the supercurrent contribution to the integral vanishes. 
The voltage is taken into account by shifting the operators Cko- by a phase exp{ieVt/h). 
Performing the integration over t' in Eq. (|T^), we find that, for small V and T = 0, we 
regain Ohm's law: J„ = R^^V. For incoherent tunneling 



8TTe^\T\^N'{0) ^ ^ t 
R„ ^ > > ^ 

while for specular tunneling. 



E E TJ^^(^±). (20a) 
q +/- ^ + ?i 



Rn' = ^El^^(q)l E (20b) 

The advantage of reporting JcRn instead of Jc is that JcRn is independent of the strength 
of the tunneling matrix element. 



IV. RESULTS AND DISCUSSION 

In this section we present the results of numerical calculations of the Josephson current 
through junction. In Fig. |^ the dependence of the critical voltage at T = on the 

chain-plane coupling is shown for an incoherent junction [Eq. (|18a|) ]. The voltage scale of 



JcRn is tens of fiV, which is two orders of magnitude lower than the value found from the 
Ambegoakar-Baratoff formulail (assuming both materials to be s-wave), and a full order 
of magnitude lower than found in the experiments of Sun et alM For small to we have a 
quadratic increase in JcRn with which is due to the distortion of the gap away from 
dx'2-y-2 symmetry by the chain-plane coupling. Increasing the coupling further, however, 
does not increase JcRn indefinitely. The maximum in the critical voltage is due to the fact 
that chain plane coupling, as well as breaking the symmetry, reduces the gap in the Cu02 
planes. This is a feature which is particular to proximity effect models. In the inset figure we 
plot the temperature dependence of JcRn for a relatively weak (to = 10 meV) chain-plane 
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coupling. The shape of the curve differs shghtly from single band models by the fact that the 
maximum value of Jc-R„ {JcRn ~ 0.026 mV) does not occur at T = 0, but at T ~ 0.4 meV. 
This happens at temperatures which are low enough that f{E^) ~ 0. As we have mentioned 
in the discussion following Eq. (0), whether is an increasing or decreasing function of 
T depends on the sign of the order parameter in the regions of the Brillouin zone which 
contribute most to the integral. Since the induced gap in the chain is smaller than the gap 
in the Pb in the region of the Brillouin zone where the chain and plane Fermi surfaces are 
far apart, the temperature dependence of JcRn at low T is determined by the induced gap. 
In Fig. |l| we can see that Aq < in the region where the induced gap is small so that JcRn 
is an increasing function of T. At larger values of T, thermal excitation of quasiparticles in 
the Pb determines the temperature dependence of the critical voltage. It is clear from this 
discussion that JcRn{T) will not have this kind of non-monotonic behaviour for an s-wave 
gap. 

In our discussion of the current through an incoherent junction, we have used the word 
"gap" in a loose sense to describe the state of the condensate. The fact that the structure 
of the gap can be distorted by the chains highlights the fundamental difference between the 
gap and the order parameter, which is defined in Eq. (^). From the definition, it is clear 
that the order parameter has the dx^-y-i symmetry of the pairing interaction regardless of 
the strength of the chain-plane coupling. In fact, our intuitive definition of the "gap" is 
more closely related to the anomalous Green's function, F, which describes both the density 
and phase of the superconducting condensate. We can rewrite Eq. (|1^) for the Josephson 
current in terms of F: 

J = ^lmY^Y.T^.nT-^.-nF'\^■,<l)FIM<l), (21) 

kq I 

where is the anomalous Green's function in the Pb, F[^ = — (Tai_k|(— ■ir)aik|(0)) is the 
anomalous Green's function in the Cu02 plane, (3 is the inverse temperature, Q = {21 + 1)tt / P 
are the fermion Matsubara frequencies and (3 is the inverse temperature. Equation ( pij ) 
makes it clear that Josephson junctions are sensitive to the structure of the condensate and 
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not the pairing interaction. In a single band material 



where E is the quasiparticle energy. In our multiband model, 

By comparing Eqs. (^) and (|2^) we can see how the chains affect the symmetry of the 
condensate, and that Fu does not share the symmetry of A^. The point we would like to 
emphasize with this discussion, then, is that a finite current through an incoherent c-axis 
junction does seem to suggest that the condensate is not (i-wave, but does not not rule out 
the possibility that the pairing interaction is (i-wave. 

In Figs. ^ (a) and (b) we plot JcRn for a specular junction, in which the tunneling matrix 
element is given by Eq. ( p.8b|) . For an isotropic Fermi surface and gap, the specular and 
incoherent cases yield identical results. As we can see in Fig. ^(a), however, the critical 
voltage is a full order of magnitude larger for a specular junction than for an incoherent 
junction. Furthermore, JcRn is a monotonically decreasing function of T. In Fig. ^(b), the 
dependence of both Jc and JcRn on to is shown. As expected, Jc vanishes as to 0, 
although here the reason is that the Fermi velocity of electrons in the 2;-direction, vl, 
vanishes. From Eq. ( pOb|) , it is clear that R~^ also vanishes as vl —>■ 0, so that the product 



JcRn is nonvanishing. This is very different from the case of incoherent tunneling where 
Rn is largely independent of to- K seems, then, that our model for specular tunneling is 
in good agreement with the observations of Sun et al^ From Fig. ^(c) we can see that for 
to = 40 meV, the critical voltage is around 1 mV, while Sun and his co-workers find critical 
voltages of 0.3-0.9 mV. This value of to is also consistent with the observed anisotropy of 
the penetration depth. AJA.,, as we have shown using a closely related model! 

The large difference between the results of the specular and incoherent tunneling is a 
reflection of the important difference between the roles of the chains in the two types of 
tunneling. As we discussed earlier, an incoherent junction is sensitive to the symmetry of 
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the condensate over the entire Brillouin zone. For a specular junction, on the other hand, 
the most important effect of the chains is to change the component of the Fermi velocity 



in favour of electrons with a large perpendicular component so that only one small region of 
the Brillouin zone, near where the chain and plane Fermi surfaces cross, contributes to the 
total tunneling current. The symmetry of the order parameter is largely irrelevant in this 
case. In the case of specular tunneling, then, the c-axis tunneling current is less a probe of 
the symmetry of the condensate than it is of the normal state band structure. 

We would like to finish this section with a brief discussion of an unresolved issue which 
is relevant to this work. Experiments on c-axis Josephson junctions have been performed 
with both twinned and untwinned crystals, and find similar values for the critical currents. 



Dynea^ has suggested that the total critical current through a junction, in which one of the 
materials is heavily twinned and has a gap which changes sign under rotations of 7r/2, should 
vanish. This is because the phase locking of the condensate at twin boundaries causes overall 
phase shifts of vr between adjacent regions. Adjacent regions should therefore have Josephson 
currents in opposite directions, and the total current should vanish. This argument is 
problematic for the work presented in this article. The argument is uncertain, however, 
because the behaviour of the order parameter at twin boundaries is not well understood. 



In this article we have calculated the Josephson current in a model c-axis YBCO/Pb 
junction. We have assumed that the YBCO is made up of alternating layers of Cu02 planes 
and CuO chains, stacked in the 2;-direction, and that the layer adjacent to the junction is a 
Cu02 plane. We take the pairing interaction in the YBCO to have d^2_y2 symmetry. For an 
incoherent junction the tunneling current is sensitive to the symmetry of the superconducting 
condensate (although not of the pairing interaction), and we find that distortions of the 
condensate due to the chains are sufficient to yield nonzero Josephson currents. The currents 



perpendicular to the junction. In Sec. |IT1| we showed that the tunneling is strongly weighted 




V. CONCLUSION 
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are an order of magnitude smaller than observed experimentally. For a specular junction, 
the Josephson current is insensitive to the symmetry of the condensate because the tunneling 
matrix element is strongly influenced by the normal state band structure. In our model, 
the Josephson current is due to one small region of the Brillouin zone over which the gap is 
roughly constant. The calculated currents for a specular junction are of the same order of 
magnitude as those found experimentally by Sun et a/.0. 
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FIGURES 

FIG. 1. Band structure of YBCO. In (a) we show the normal state Fermi surface for the 
dispersion given in Eq. (p!l|). The Fermi surface has two different pieces since the unit cell consists 
of a chain and a plane. The Fermi surfaces are shown for a range of chain-plane coupling values 
(0 < t < 50 meV). In the absence of chain-plane coupling (i = 0), the Fermi surfaces cross. As 
t is increased, the Fermi surfaces are pushed apart. The dashed line is the line along which 



vanishes. In (b) and (c)we plot the quasiparticle excitation energies E± [Eq. (8b)] along ky = and 
ky = 8kx/3 respectively. We have taken t = 40 meV for these curves and have plotted t = limits 
of E± [(^f + A^)-*^/^ and |^^.|] for comparison. At temperatures lower than the of Pb, thermal 
excitation of quasiparticles in the YBCO is limited to the nodes of E^. All results presented in 
this paper are for cji = 100 meV, (T2 = 60 meV, /ii = —80 meV, /i2 = 40 meV. 

FIG. 2. Critical voltage, JcRn, for an incoherent junction. We show the dependence of 
JcRn{T = 0) on the chain-plane coupling parameter to- For small to we find the expected quadratic 
increase in the critical voltage as the coupling distorts the gap away from (1^2 _y2. At larger to, 
the coupling to the chains weakens the condensate, as well as distorting its structure, leading to 
a maximum in the curve. This curve shows that, for the simple model of incoherent tunneling, 
distortion of the condensate by the chains cannot account for the experimentally measured values 
of JcRn ^ 0(0.5 mV). The inset shows the temperature dependence of the critical voltage for 
to = 10 meV. The curve increases slightly at low T. This stems from the multiband nature of the 
YBCO and the antisymmetry of the order parameter. 
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FIG. 3. Critical voltage, JcRn, for specular tunneling. Here the tunneling matrix element is 
chosen to conserve the component of the wavevector parallel to the junction. Furthermore, the 
tunneling matrix element is weighted in favour of particles with a large perpendicular velocity [Eq. 
(|18b| )]. In (a), we show the temperature dependence of JcRn for our multiband d-wave model 
(solid line), and for the simple s-wave model of Ambegoakar and Baratoffii(dashed line). In order 
to make the T = values of JcRn agree, we have taken 2A/Tc ~ 0.12 in the s-wave model. The 
most important difference between this figure and Fig. ^ is that critical current is a full order of 
magnitude larger here. In (b), both JcRn{T = 0) and Jc(T = 0) are plotted as functions of to- 
The magnitude of Jc{T = 0) is arbitrary since the tunneling probability, P, has not been specified. 
The critical voltage is nonvanishing as to ^ because Rn diverges. 
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